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We report on a search for the decay Kl — > it°fJ- + fl~ carried out as a part of the KTeV experiment 
at Fermilab. This decay is expected to have a significant CP violating contribution and a direct 
measurement will either support the CKM mechanism for CP violation or point to new physics. 
Two events were observed in the 1997 data with an expected background of 0.87 ± 0.15 events, and 
we set an upper limit B(Kl — > 7r°/i + /i~) < 3.8 x 10~ 10 at the 90% confidence level. 

13.20.Eb, 11.30. Er, 14.40.A 

The decays Kl — > 7r ^ are interesting decays for the study of CP violation and can be used to search for new 
physics. There are three expected contributions to the amplitude: a CP conserving contribution which proceeds 
through the 7r°7*7* intermediate state, an indirectly CP violating contribution from K\ — ► tt°U, and a directly CP 
violating contribution from electroweak penguin and W box diagrams jl] -0]. Branching ratio predictions in theories 
containing exotic (e.g., SUSY) particles that contribute to the penguin amplitudes are significantly higher Q. 

The sizes of the three contributions depend on the flavor of the final state lepton. The greatest theoretical interest 
is in the K L — > tt°vv case, where the direct CP violating amplitude dominates and a theoretically clean measurement 
H of the Wolfenstein || parameter r] should be possible. However, measuring a final state with two neutrinos and 
a single pion is experimentally challenging, and the current experimental limit remains four orders of magnitude 
above the Standard Model expectation of ~ 3 x 10~ n . 

In contrast, Kl — > 7r°e + e~ and Kl — ► 7r°/i + /i - are comparatively straightforward to detect, although all three 
amplitudes are present in these modes. This letter presents a new limit on B(Kl ~ * 7r°/j + /i~ ); the existing limit || 
is 5.1 x 10~ 9 at the 90% C.L. For the data taken by KTeV in 1997 and discussed here, a single event observed in the 
muon mode would correspond to a branching ratio of 7 x 10 -11 , which approaches the Standard Model expectation 
|,[cj of B(K L -> ttV+zO ~ (0-44 - 1.00) x 10~ n . 

Figure [l] shows a plan view of the KTeV detector, which has been described elsewhere [0,0- An 800 GeV proton 
beam, with typically 3.5 x 10 12 protons per 19 s Fermilab Tevatron spill every minute, was targeted at a vertical angle of 
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4.8mrad on a 1.1 interaction length (30 cm) BeO target. Photons were converted by 76 mm of lead immediately down- 
stream of the target. Charged particles were then removed with magnetic sweeping. Collimators defined two 0.25 /isr 
beams that entered the KTeV apparatus 94 m downstream of the target. The 65 m vacuum (~ 10~ 6 Torr) decay region 
extended to the first drift chamber. The spectrometer consisted of a dipole magnet surrounded by four (1.28 x 1.28 m 2 
to 1.77 x 1.77m 2 ) drift chambers with ^100 /.tm position resolution in both horizontal and vertical directions. Helium 
filled bags occupied the spaces between the drift chambers; the magnetic field imparted a ±205MeV/c horizontal mo- 
mentum kick. The spectrometer had a momentum resolution of a(P)/P — 0.38% 0.016%P, where P is in GeV/c. 
The electromagnetic calorimeter consisted of 3100 pure Csl crystals. Each crystal was 50 cm (27 radiation lengths, 1.4 
interaction lengths) long. Crystals in the central 1.2 x 1.2 m 2 section of the calorimeter had a cross-sectional area of 
2.5 x 2.5cm 2 , and those in the outer region (out to 1.9 x 1.9m 2 ) had a 5 x 5 cm 2 area. The calorimeter's energy reso- 
lution for photons was a(E)/E — 0.45% © 2%/^/E, where E is in' GeV, and its position resolution was ~1 mm. Nine 
photon veto assemblies (lead scintillator sandwiches) detected particles leaving the fiducial volume. Two scintillator 
hodoscopes in front of the calorimeter were used to trigger on charged particles. The hodoscopes and the calorimeter 
had two holes (15 x 15 cm at the calorimeter) to let the neutral beams pass through without interaction. The muon 
filter, located behind the calorimeter, was constructed of a 10 cm thick lead wall followed by three steel walls totaling 
511cm thickness. Scintillator planes with 15 cm segmentation in both horizontal and vertical directions (MU3) were 
located after the third steel wall. The segmentation was comparable to the multiple scattering angle of 10 GeV muons 
at MU3. Pion punch-through probabilities, including decays downstream of the calorimeter, were taken as a function 
of momentum from Kl — > 7r ± e =F ^ data and are on the order of 2 x 10~ 3 . The data acquisition system reconstructed 
events online, and the results were used to filter the data. 

The signature we searched for is two tracks from oppositely charged particles with a common vertex that deposit 
little energy in the calorimeter and created two hits at MU3 from these muons. The 7r° creates two electromagnetic 
showers in the calorimeter with m 77 = m,o and which are unassociated to tracks. 

There are three important backgrounds. The first is Kl — > 7r + 7r~7r where both either decay upstream of the 
calorimeter (decay in flight) or punch through to MU3. The second is Kl — > n^/x^u with one decay in flight or punch- 
through and accidentally coincident calorimeter activity that appears as a tt°. The third and largest background is 
the radiative muonic Dalitz decay Kl — ► M + / i_ 77 when m 77 = m OT o. Because of the low expected ]l2] ] branching 
ratio, Kl — * 7r°7r ± /i :F i/ is not a large background. 

Two triggers were used for this analysis. To determine the number of Kl decays in the data, we identified 
Kl — > 7T + 7r _ 7r° decays in a minimum bias trigger. This trigger required hits in the trigger hodoscopes and the drift 
chambers which were consistent with two coincident charged particles passing through the detector. Events with a 
reconstructed vertex from oppositely charged tracks were recorded with a prescale factor of 500:1. For the signal 
trigger further requirements were made. Two or more hits in MU3 were required, and activity in the photon veto 
counters rejected events. The trigger system counted the number of calorimeter clusters over ~ 1 GeV in a narrow 
(20nsec) time gate; for the signal mode, at least one such cluster was required. We also required that the calorimeter 
energy reconstructed online and associated with each track be less than 5 GeV. The signal trigger was not prescaled. 

Muons passing through the calorimeter typically deposit ~400MeV, and so in searching offline for the signal we 
required that the clusters associated to the tracks had less than 1 GeV of energy. Track momenta were required to be 
greater than 10 GeV/c to ensure that they penetrated to MU3 and less than 100 GeV/c to ensure that the momentum 
was well measured. The two-muon system was required to have a mass less than 350MeV/c 2 to reduce backgrounds 
from Kl — ► n ± /i T u. We required two non-adjacent hits in both views of MU3. We did not compare the MU3 hit 
positions to the extrapolation of the tracks to MU3, because the major backgrounds passed this requirement as well 
as the signal did. 

To suppress ir^ decays in flight, we required that the reconstructed vertex occurred in the beam volume of the 
decay region and had a x 2 of 10 for 1 d.o.f. or less, and that the track segments upstream and downstream of the 
spectrometer magnet passed within 1mm (94% signal acceptance) of each other at the bend plane of the magnet. 
The mass of the two unassociated clusters under the hypothesis that they were produced by photons from the decay 
vertex was required to be between 135±6MeV/c 2 (±2.5<r). These clusters both had to have been found by the trigger 
cluster counter. 

A number of kinematic criteria were studied to suppress background from Kl — > /i + /i~77 decays. The kinematics 
of this decay are very different from the analogous Kl — > e + e~77 background to Kl — ► ir°e + e~ , and the branching 
ratio is lower. Consequently the methods |l^,|l4| which are effective in the case are less effective in the fi^ case 
and were not applied. 

It is possible to suppress backgrounds from Kl —* 7r + 7r~7r°, and Kl - > n^fi^v by using 
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where mx is the kaon mass, m w the mass of the two-muon system, m^a the ir° mass, and is the two-muon 

systems's momentum perpendicular to the kaon momentum. This quantity is proportional to the 7r° momentum 
squared in the Kl flight direction in the frame colinear with the Kl but where the pair has no longitudinal 

momentum. We required Rj^ to lie between -0.01 and 0.10 GeV 2 /c A . This cut keeps 89.2% of the signal and rejects 
73% of Kl — > 7r ± ^i T ^ decays with coincident photons and 95% of the Kl — > 7r + 7r~7r decays. 

To ensure that we observed all the products of a Kl decay, we required that the total squared momentum transverse 
to the Kl flight direction (P]_) be less than 100 MeV 2 /c 2 , and that the reconstructed mass (m) of the Kl be between 
492 and 504MeV/c 2 . With these requirements, which were selected by examining Monte Carlo simulation results and 
data outside the signal region before examining the data for Kl — * 7r°/i + /i~ candidates, the overall acceptance for the 
signal was 5.0%. The simulation distributed the products of the decay uniformly in phase space. 

Figure || shows the P± 2 vs. m distributions for the data, and Fig. || shows the mass distribution after the Pj_ 
requirement for the data and the backgrounds as estimated from the simulation. The correspondence between the 
data and the simulation is good, and is also good in the distributions of m MAl , Pj, track momentum, vertex position 
and (for Kl — > 7r + 7r~7r°) of m 11 . The background levels in the signal region are given in Table |[ they are calculated 
from the simulations, published |16[ branching ratios, and the number of Kl decays in the data sample. Although 
we have observed Kl — > /x + /i~77 ]15| in this data, this background was more precisely estimated from QED and the 
measured B(Kl — ► [i + (J,~-f); a value of B(Kl — > n + fJ.~"fy) = (9.1 ± 0.8) x 10~ 9 with an m 77 > 1 MeV/c 2 cutoff was 
used. All Monte Carlo samples were over 10 times the data sample. The background from Kl — > 7r + 7r~ + 27 (Acc) 
and Kl — * 7r + 7r~7 + 7 (Acc) was negligible. 

To normalize any possible signal's branching ratio, we identified Kl — > 7r + 7r~7r decays in as similar a manner to the 
Kl — > 7r°/i + /i~ identification as possible. Apart from the trigger differences, the calorimeter energy requirement for 
clusters associated to tracks was changed to less than 0.9 times the momentum measured with the spectrometer; the 
MU3 and Rj^ requirements were removed, and m w ± rather than m^± was used in calculating kinematic quantities. 

The acceptance for the normalization mode was 8.1%. There were (268 ± 0.4stat ± 0.4mc ± 4.3br) X 10 9 Kl decays 
between 90 and 160 m from the target with Kl momentum between 20 and 220GeV/c. 

In calculating the number of Kl decays in the data and the acceptance for Kl — * 7r /! + /i~ relative to the acceptance 
for Kl — > 7r + 7T~7r°, we allowed for the uncertainties summarized in Table [n]. We calculated the Kl flux using 
Kl — > 7r ± /i T ^ rather than Kl — > tt + tt^tt° decays and attributed the difference of 4.20% to the quality of our simulation 
of muons in the detector. We varied the scale and resolution of the calorimeter and spectrometer in the simulation to 
conservatively cover the range of variations seen in the data. Apart from uncertainties in published branching ratios, 
other sources of uncertainty were small. 

From Figs. || and ||, two events exist in the signal region for the data. Sidebands in both m and P± 2 show 
correspondence between data and background predictions as given in Table III , With the above acceptance and Kl 
flux, and allowing for a background level of 0.87±0.15 events, we set an upper limit B(Kl — > 7r°/i + ^~) < 3.8x 10~ 10 
at the 90% confidence level. 

This limit is approximately one order of magnitude more stringent than the previous limit. If we assume that 
the only contribution to the branching ratio is from direct CP violation, we may conclude that |»;| < 7 at the 90% 
confidence level. In comparison to the search for Kl — > 7r°e + e~, Kl — + 7r°/i + /i~ searches will have better single event 
sensitivity for any given sample of Kl decays because the level of irreducible Kl — * ? + /~77 background is less. While 
not yet as sensitive as B decays, where a recent indirect global analysis |ll| finds 77 to be below 1, it is valuable to 
test if the same parameterization is valid for both B and K decays. 
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TABLE I. Summary of expected backgrounds in the sig- 
nal region. "Acc." refers to particles from other sources 
which were accidentally time-coincident with the interesting 
decay; "D" and "P" refer respectively to decay in flight or 
punch-through. Limits are 90% C.L.; uncertainties are due to 
uncertainties in published branching ratios, simulation statis- 
tics, and the statistics of the normalization mode. 





Decay mode 




Expected number of events 


K L - 






0.373 ± 0.032 


K L - 


-> ^"V~7 + 7(Acc) 




< 0.029 


K L - 


-> TT + TT-TT (DD) 




0.252 ± 0.095 


Kl - 


■* TT + TT-TT (DP) 




0.007 ± 0.007 


K L - 


-» TT + TT'TT (PP) 




0.007 ±0.007 


K L - 


-» Tr+zi^ + 27(Acc) 


(D) 


0.161 ±0.093 


Kl - 


-> ■K ± [i T v + 27(Acc) 


(P) 


0.063 ± 0.037 


K L - 


-> tt'V+^T^ (D) 




0.009 ± 0.009 


K L - 


-> TT°TT+/i+i / (P) 




< 0.009 




Total 




0.87 ±0.15 
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TABLE II. Systematic and statistical sources of uncer- 
tainty. Sources marked with (*) contribute to uncertainty 
in both the Kl flux and the acceptance for Kl — > tv°h + /j,~ 
relative to the acceptance for Kl — » 7r + 7r~7r°; other sources 
contribute only to the acceptance ratio. 



Source Relative Uncertainty 



B(K L TT+TT-TT ) 


1.59% (*) 


Data statistics for K L — > 7r + 7r~7r () 


0.16% (*) 


Simulation statistics for Kl — > 7r + 7r~7r° 


0.14% (*) 


Simulation statistics for Kl — > ir n + fi~ 


0.16% 


Calorimeter scale and resolution 


3.33% 


Spectrometer scale and resolution 


1.12% 


Muon identification 


4.20% 


Signal trigger requirements 


0.80% 


Vertex quality requirement 


0.22% 


Spectrometer wire inefficiency 


0.15% 



Total ' 5.77% 



TABLE III. Number of observed and predicted events in 
regions near the signal region. 

Region Prediction Observed 

480MeV/c : '< m <492MeV/c ;i ; PJ<0.0001( GeV/cf 2.75±0.73 2 

504MeV/c 2 < m <515MeV/c 2 ; <0.0001( GeV/c) 2 0.14±0.08 

492MeV/c 2 < m <504MeV/c 2 ; 0.0001( GeV/c) 2 <P|<0.0005( GeV/c) 2 2.65±1.53 2 

492MeV/c 2 < m <504MeV/c 2 ; 0.0005( GeV/c) 2 <P|<0.0010( GeV/c) 2 2.80±2.59 2 
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FIG. 1. The KTeV detector configuration for rare decay 
studies. The TRD, Hadron Veto, and the Photon Veto in the 
beamline were not used in this analysis. 
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FIG. 2. Reconstructed P ± vs. m after all other selection cri- 
teria, for the data. The box indicates the signal region; events 
in the lower left are predominantly Kl —* tv + tv~tv° with tt^ 
decays in flight. 




FIG. 3. Reconstructed m after all other selection crite- 
ria, for the data (dots), and the total background estimate 
(line). The arrows indicate the range of values of m accepted. 
The small bump in the background at the Kl mass is from 
Kl — * /i + /i~77. 



G 




0.46 0.48 0.5 0.52 0.54 



m (GeV/c 2 ) 



